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Abstract 

Orbital Angular Momentum (OAM) twist-two parton distributions are studied. At the 
low energy, hadronic, scale we calculate them for the relativistic MIT bag model and for 
non-relativistic potential quark models. We reach the scale of the data by leading order 
evolution using the OPE and perturbative QCD. We confirm that the contribution of 
quarks and gluons OAM to the nucleon spin grows with Q 2 , and it can be relevant at 
the experimental scale, even if it is negligible at the hadronic scale, irrespective of the 
model used. The sign and shape of the quark OAM distribution at high Q 2 may depend 
strongly on the relative size of the OAM and spin distributions at the hadronic scale. 
Sizeable quark OAM distributions at the hadronic scale, as proposed by several authors, 
can produce the dominant contribution to the nucleon spin at high Q 2 . As expected by 
general arguments, we obtain, that the large gluon OAM contribution is almost cancelled 
by the gluon spin contribution. 
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1 Introduction 



Understanding how the partons carry the angular momentum in the nucleon has be- 
come a main effort of present day physics. The quark spin contribution AE is well defined 
in QCD [||], thus measurable, and the gluon spin contribution Ag is measurable, despite 
some problems of definition || and therefore both have been studied intensively. The 
situation with respect to the quark, L q and gluon, L g , Orbital Angular Momentum (OAM) 
is less satisfactory. 

It is well known, that the most natural definition of OAM for quarks and gluons cannot 
be separated in a gauge invariant way from the corresponding spin terms @. However, 
recently, new definitions of angular momentum have been implemented to accommodate 
gauge invariance j|, ||, and from them adequate twist two OAM distributions have been 
constructed. 

In this respect, three different possibilities have been investigated. One proceeds by 
choosing a particular gauge ||, |[ 0. The OAM operator leads to the naive definition of 
orbital angular momentum up to effects not controlled by the gauge fixing. A second, 
maintains gauge invariance, by loosing covariance, defining the distributions only in the 
class of reference frame where the nucleon has a definite polarization [Q. In this case 
the resulting distributions can be related to the forward limit of off-forward quantities, 
and are measurable [|, ^ [| . The last proceeds by defining OAM operators such that the 
distributions are gauge invariant [§] . Furthermore, in the light-cone gauge, they reduce 
to the natural definitions |§. At present, however, no physical process has been proposed 
to access them. 

Evolution equations for the OAM distributions have been derived || [7], [8], [U], , and 
it has been shown || that: 1) at leading order (LO) the evolution is determined by known 
(polarized and unpolarized ) anomalous dimensions ; 2) at LO the evolution equations for 
the different definitions of the operators coincide; 3) for the second procedure above, the 
evolution is related to that of the polarized and unpolarized distributions at any order. 

The evolution equations have been numerically solved by using as input data-inspired 



OAM distributions [12 



During the last few years we have developed a scheme to study distributions based on 



model calculations [fL3| . This procedure has been developed thus far only to leading order 
in the twist expansion and therefore we should compare with the data only at high Q 2 , 
where the contribution from the non leading twists vanishes. 

We have analyzed the quark and gluon spin distributions, as well as odd chirality 
ones [H|, in this scheme. In this Letter we present the OAM distributions obtained by 
evolving those, properly calculated within relativistic and non relativistic models, to the 
experimental scale. 

We present in the next section the theoretical framework behind our calculation and 
discuss our definitions of the OAM distribution and their relation to the model distribu- 
tions. The following section deals with the results of the evolution for several scenarios at 
the original hadronic scale. We analyze how hadronic structure influences the DIS results. 



1 



Finally we extract in the conclusions the phenomenological implications of our study. 



2 The theoretical framework 

It has been suggested in the past that the OAM contribution to the nucleon spin can be 
large at low [16| , as well as large energy scales. Since the quark spin, A£(Q 2 ), 
and the gluon spin, Ag(Q 2 ), are observables, one may use the Spin Sum Rule |§ for the 
nucleon, 

^A£(Q 2 ) + Ag(Q 2 ) + L q {Q 2 ) + L g {Q 2 ) = \ (1) 

to determine the global OAM contribution to the spin. The problem of separating the 
OAM in the quark and gluon fractions in a gauge invariant way, already addressed in || , 
is a cumbersome one. In principle it has been solved in ||, where a proper definition of 
the distribution functions for every term in ([!]) has been introduced. Unfortunately these 
distributions have not yet been related to observable quantities. Another approach to it 
has been proposed in Ref. ||, where the authors have shown that a gauge independent 
decomposition for the total angular momentum, J, in quark, J q , and gluon, J g , contri- 
butions exists in an interacting field theory. Moreover, J q and J g are gauge independent, 
and therefore they can be measured, for example, through a forward limit of off-forward 
quantities []. Based on these developments, a gauge independent, twist-2, new definition 
for the quark OAM distribution, whose evolution equation is known, has been given in 
ref. |J. This definition mixes the polarized and unpolarized singlet quark distributions 
with the forward limit of OFFPD. Since the OAM matrix element is not a Lorentz scalar, 
there is an ambiguity in this definition for any relativistic quantum theory. To avoid it, 
one must take for QCD a system of coordinates where the nucleon has a definite helicity, 
thus loosing covariance. 

We proceed to use this last theoretical development, to perform a phenomenological 
model analysis of the OAM distributions, using relativistic as well as non-relativistic quark 
models. 

Let us first discuss the definition of the OAM for a relativistic model. In this case we 
follow the definition provided by ref. j8|, which is valid, in a class of coordinates in which 
the nucleon has a definite helicity. If we assume, for example, that the nucleon is moving 
in the z direction and is polarized with helicity +1/2 , the quark OAM distribution is 
given by, 

L q (x, Q 2 ) = \ Q 2 ) + E q (x, Q 2 )) - A£(x, Q 2 )] , (2) 

where Q 2 ) (A£(x, Q 2 )) is the usual unpolarized (polarized) singlet quark distribution 
and E q (x,Q 2 ) is the forward limit of the helicity-flip, chiral odd, twist-two OFFPD 

1 A thorough discussion of off-forward parton distributions (OFFPD) can be found in Refs. jl?], Q. 
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E(x, A 2 , A • n) 0, 17 1 . The latter quantity is defined through the twist- 2 part in the twist 
expansion for the following light cone correlation function, 



^<P^(-An/2) 7 ^(An/2)|P> = 

H(x,A 2 ,A-n)U(P')YU(P) + E{x, A 2 , A • n)U(P 



2M 



U{P) 



(3) 



where P(P') is the 4- momentum of the initial (final) nucleon in a virtual Compton scat- 
tering process, A /t = P' — P^, n = (1, 0, 0, — 1)/2A, A is fixed by the choice of the reference 
frame, and H(x, A 2 , A • n) is the helicity-conserving, chiral even, twist-2 OFFPD whose 
forward limit is the usual forward unpolarized parton distribution. In the definitions 
above, the forward limit corresponds to A 2 — > 0, A • n — > 0. 

In the non-relativistic case, the nucleon wave function is given in general by an ex- 
pansion in terms of the eigenstates of some approximate hamiltonian. Let N label the 
quantum numbers of the eigenstates, i.e., principal quantum number, orbital angular 
momentum, spin, then 



V(Pl,P2,P3) = J2 a N^N(PuP2,P3) 
N 



(4) 



The OAM quark parton distribution, generalizing our approach developed for unpolarized 
and polarized distributions |T3]|, is determined by, 



where 



roo 

L q (x,Q 2 ) = 2irMY,M 2 / d PP L z qN (p) 



L q,N(p) = (^N(Pl,P2,Ps)\J2 L i S (P~Pi)\^N(Pl,P2,P3)}, 

i=l 



and the lower integration limit is 



\p-i x )\ = y 



X 



m q \ 2 1 
~M) x 



(5) 



(6) 



(7) 



Having set down the framework which defines the OAM distributions at the hadronic 
scale, we proceed to calculate them explicitly in two models and to study their evolution. 



3 Analysis of Orbital Angular Momentum distribu- 
tions at high Q 2 

We proceed to study the OAM distribution in two different scenarios for proton structure: 



i) a non-relativistic scheme based on the Isgur-Karl model |19| ; ii) a relativistic approach, 



as described by the MIT bag model [20 
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For the non-relativistic scenario of proton structure we consider initially the Isgur-Karl 
model with a proton wave function given by a harmonic oscillator potential including 



contributions up to the 2hu shell In this case the wave function, Eq. (|J), is given 
by the following admixture of states 

\N) = a s \ 2 S 1 / 2 )s + 0'S'\ 2 S' 1 / 2 )s + o,m\ 2 Si/2)m + «d| 4 -Di/2)m , (8) 

where we have used the spectroscopic notation \ 2S+1 Xj) t , with t = A,M,S being the 
symmetry type. The coefficients were determined by spectroscopic properties to be: a$ = 
0.931, a s > = -0.274, a M = -0.233, a v = -0.067. 

Calculating Eq. @ using the wave function Eq. one gets for the parton distribu- 
tion, Eq. (|), 

T t 2x | |2 M /3\ 3 / 2 /lpi 13p 2 23\ !4 

where a 2 is a parameter of the model, and p-(x) is given by (|7|). Note that only the small 
| 4 -Di/2)m wave component gives a contribution to this OAM distribution. 

In the relativistic scenario, we proceed to evaluate in the MIT bag model the twist-two 
distribution Eq. @. The term AS (a;,// 2 ,) in Eq. (0) is discussed for the bag in ref. f2~4} . 

To calculate Eq.(|2|) we need the quantity S(x, (/%) + E q (x,fj, 2 ) ), which is given by the 
forward limit, A 2 -> 0, A-n -> 0, of the OFFPD distribution H(x, A 2 ,A-n) + E(x, A 2 , A- 
n). The result for the latter in the bag, calculated by using the general definition (Eq. 



)), is to be found in Eq.(29) of ref. [p2 |, from which we have obtained the forward limit, 



Z(x,nl)+E q (x,n 2 ) = Z 2 (0) / dk ± d(f>k ± {tl(k)+ (10) 



+ MMx-i)(2,-x + i)^} 



where 



iV^ = ^ (11) 

2i?3(^ - l)j'o 2 M 1 J 

Z(0) = N 2 J q R drr 2 ( Jo 2 (^r) + j{ (^r)) (12) 

u; = 2.04 is the lowest energy eigenvalue for the quark in the bag, j (i) are the spherical 

Bessel functions, k is the momentum of the quark in the bag, k' = \k\, k' = k + 77 A, rj is 
a parameter of the calculation in |22| and, eventually 

= j (q;o) cos(kR) - j (fc.fi) cos^p 

° l J u 2 -pR? 1 J 



and 
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and 

1[ ' ~ u 2 -k 2 R 2 ' 1 ' 

We recall, that the calculation of H(x, A 2 , A-n)+E(x, A 2 , A-n) in [22], whose forward 
limit we are using here, is performed in the Breit reference frame. The bag model has 
no translational invariance and therefore the choice of frame must be associated to the 
model assumptions. Our scheme requires that the evaluation is performed in the rest 
frame, thus we have to fix here the parameter rj to the value of 0.75 which corresponds, 



according to Ref. p2| , to the unboosted calculation. 

In both cases, Isgur-Karl (IK) and MIT, we use the corresponding support correction 
as defined in |y| and [p3[j , respectively. 

Once the distributions L q (x,/j,l) have been obtained at the low (hadronic) scale of 
the model, we perform a LO QCD evolution according to the equations displayed in refs. 
0, H |l(| [Tl| and numerically solved in \12\. These contain a complicate mixing between 



L q ( g )(x, ^q), AS(x, and Ag(x, /Xq), a feature which will be very relevant in the analysis 
of the data. 

The results of our analysis are shown in Figs. 1 through 5. 

In Fig. 1 we show the IK result for quarks (gluons) in (a) ((b)); The full curve 
corresponds to the initial distribution, which is missing in (b), since we start from a very 
low hadronic scale where no valence gluons exist. At LO the hadronic scale corresponds 
to fiQ ~ 0.08 GeV 2 . The dashed curves correspond to the result of the evolution from the 
hadronic OAM distributions, to Q 2 = 10 GeV 2 (short-dashed) and to Q 2 = 1000 GeV 2 
(long-dashed). 

We can summarize the results of the calculation as follows: i) the evolved distributions 
are negative; ii) the magnitude of the distributions increases with Q 2 at low x; iii) the 
magnitude of L q (x, Q 2 ) is small but increases with respect to the tiny starting distribution 
(Eq. (I)); iv) though gluons are assumed to be negligible at the hadronic scale, L g (x, Q 2 ) 
becomes much larger than L q (x, Q 2 ) at high Q 2 . 

We show in Fig. 2 the same analysis for the MIT bag model. Here the initial L q is 
much larger, but the result of the evolution is qualitatively basically the same. 

Thus our first conclusion is that there is little model dependence for different initial 
OAM distributions. In this case their structure does not seem to influence very much the 
evolution. It is clear that the other inputs of the equations, AE(x, /Zq) and Ag(x, /Iq), are 
the dominating features. 

For these results the initial scale has been very low and therefore no initial gluon 
distribution was required. In Fig. 3 we show the result of the evolution for the IK model 
starting from a higher scale, /1q ~ 0.23 GeV 2 . At this hadronic scale, for LO, about 40 % 
of the proton momentum must be carried by the gluons. The polarized gluon contribution 
is built starting from the valence quark distributions as done in (TBj and suggested in J25J, 
and we define L g (x,fil) from L q (x,fil) using the same prescription. 

Again the same features as in the first analysis are found. The initial gluon distribution 
is so small, in relation with the final one, that it does not show up in (b). 
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We note that this choice of the initial scale is the same adopted in the numerical 



analysis of ref. |T2|]. An important check of the consistency in our evolution is that we 
have recovered the results of this reference by inserting their initial distributions in our 
code. 

Recapitulating, our results show that the input distributions L q (x, //q) and L g (x, do 
not seem to determine the behavior of their evolved ones, which turns out to be governed 
by the the singlet polarized distributions Ag(x, /1q) and AS(x, /ig), due to their mixing 
in the evolution equations. To check to what extent such a statement, stressed also in 



12], is valid, we analyze in Fig. 4 the results of a modification of the IK model, the so 



called D model, already studied in ||26|| . In this variant model, the D wave probability 



is set large to reproduce the axial coupling constant of the nucleon |I5| . This condition 
requires the following choice of the parameters in Eq. (H): as — 0.894, as' = 0, = 0, 
a-p = —0.447, i.e., the probability to find a nucleon in the D wave is about 20 %. Moreover 
we have introduced also in the D Model scenario polarized valence gluons, as we did for 
the IK scenario of Fig. 3. From the figure it is clear that, while the result for the gluons 
does not differ in a relevant way from the ones found before with the various models, the 
result for L q (x, Q 2 ) does, this distribution becoming rather large and positive for large x. 
It is important to stress that in order for this to occur two mechanisms were needed, a 
large initial OAM distribution (as in the MIT bag Model), and a higher hadronic scale 
(as provided by the valence gluons). As shown previously, the independent action of 
the two mechanisms does not lead to this behavior. It is the joint action of large initial 
OAM distributions and a softer evolution which is responsible for it. Thus, Ag(x, /Zq) and 
A£(x, iIq) are governing the evolution as long as they are much larger than L q (x, /z 2 ,) at 
the initial scale. When they have similar size the above statement is not true any more. 
Note that the IK interaction together with the choice of parameters of the D model does 
not describe the hadron spectrum. Nonetheless, other models of interaction (for example 



T5|, [L6|]) predict at the low scale a large OAM and fit the spectrum. We conclude therefore 
that a precise knowledge of the OAM distributions will serve to distinguish among the 
models. 

In Fig. 5 we show the evolution of the various contributions to the spin sum rule, 
Eq. (JI|). We write here below the LO-evolution equations for the first moments of the 
distributions involved, corresponding to three active flavors (cf. 



&9(Q 2 ) = i (b- 1 - l) AE(tfj) + b^Ag^l) , (16) 



AS(Q 2 ) = AE(/Xq) , (15) 
4 
9 

L q (Q 2 ) = (6-S-l)iAEM+6-fiL 8 0ig)-|(6-e - l) , (17) 
L g (Q 2 ) = 6-tA^/i 2 ) - Ag(Q 2 ) + b^L^) - A ( r f _ i) , (1 8 ) 

where b = a s {Q 2 )/a s (^l). 

Fig. 5 (a) corresponds to the modified IK scenario used in Fig. 3, where at the scale of 
the model the OAM carried by quarks and gluons is very small and the rest of the proton 
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spin is almost equally shared between quarks and gluons spins (/^(/Iq) + Lg(i$) ~ 0.01, 
AS(/!q) ~ 0.48 and Ag(^) ~ 0.25), whereas Fig. 5 (b) corresponds to the extreme 
scenario used already in Fig 4, the so called D model, with a large initial OAM (L 9 (/Zq) = 
0.145, Lg(fil) = 0.055, AS(^o) ~ 0.4 and Ag(^l) ~ 0.1). As predicted by total angular 
momentum conservation |l| and already obtained in |12] as a model-independent feature 
of the evolution equation , it is seen that at large Q 2 the huge negative contribution 
L g (Q 2 ) basically cancels out with the positive Ag(Q 2 ). Anyway, the role of the quark 
OAM is found to be very important in the second scenario (cf. Fig. 5(b)), being at large 
Q 2 the largest contribution to the saturation of the spin sum rule. Again, we see that 
quark OAM, due to evolution, can be important at large Q 2 if it is not negligible at the 
scale of the model, whereas the gluons OAM, though it is large, is basically cancelled by 
the gluons spin. 

Before concluding this section a comment on our perturbative evolution is advisable. 
As has been discussed previously the scale of the models is rather low. It is therefore 
natural to think that NLO corrections are relevant and could change the interpretation 
of the results. Actually, our past experience [|13|, [14| in similar situations suggests that 
NLO corrections are at the level 30 — 40% at most, and moreover, they never spoil the 
LO behavior. We are thus confident that the main results of our analysis will survive, at 
least in a qualitative fashion, after the implementation of NLO corrections. 



4 Conclusions 

We have studied the OAM distributions as defined newly in order to take into account 
gauge invariance. Our study, as in previous occasions fl3f , is at the heart of the ill-labelled 



spin crisis. We define a hadronic scale at which quark models should be operative. We 
use the quark model calculations as non-perturbative input for the evolution equations 
and therefore to predict the experimental outcome at the DIS scale. We have seen in this 
way that much of the folklore, associated with the spin crisis, has banished. However the 
procedure has served to distinguish, in detailed calculations, among the different models. 
In here we perform an analysis of the OAM distributions. 

We have seen that evolution, as in the other calculations, plays a major role in the 
outcome of the predictions. The fact that the gluon and quark spin singlet distributions 
mix in these equations with the quark OAM distributions, implies that for large Q 2 large 
contributions from the OAM are to be expected, even if they are not present at the 
hadronic scale. Thus, two scenarios arise in a natural way. One, the more conventional 
one, as described by the more traditional models, is defined by quark OAM distributions 
at a small hadronic scale. In this scenario the evolved distributions are large and negative 
and almost model independent and angular momentum DIS physics is dominated by the 
quark and gluon spin singlet distributions, not by OAM distributions at the hadronic 
scale. A second scenario is defined by quark and gluon OAM distributions at a higher 
hadronic scale. In the latter, soft evolution scenario, the initial distributions are important 
and therefore DIS physics may be able to discriminate between different models. If the 
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OAM distributions are large the outcome of the evolution is very strongly dependent on 
the initial distributions and completely different from that of the first scenario. 

Finally the gluon OAM distributions become huge through evolution, even if they are 
not present at the hadronic scale. However, as it is well known |I| , the gluon OAM and 
gluon spin contributions cancel to a great extent in the nucleon spin, but not so in other 
moments. 

Our past experience suggests that LO results provide a reasonable qualitative approx- 
imation and we do not expect that NLO corrections can spoil their general features. 

It is clear that many phenomenological implications have arisen of our study. A careful 
analysis of gauge invariance [|], ||, || has permitted us to obtain many observables, which 
may not only lead to a better understanding of the proton spin, but to describe the 
proper behavior of QCD at low energies, i.e., in the confining region. These observations, 
implemented in our scheme, are instrumental in defining the picture of the proton at the 
hadronic scale that should be used for describing low energy properties. 
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Captions 



Figure 1: Proton OAM distributions in the IK model for the quarks, (a), and for the glu- 
ons, (b). The full curves show the initial distributions at the hadronic scale of n 2 = 0.08 
GeV 2 , where a negligible fraction of the nucleon momentum is carried by the gluons; the 
dashed curves represent the LO-evolved distributions at Q 2 = 10 GeV 2 ; the long-dashed 
curves give the LO-evolved distributions at Q 2 = 1000 GeV 2 . 

Figure 2: The same as in Fig.l, but for the MIT model. It should be noticed that 
the distribution at the hadronic scale of = 0.08 GeV 2 is shown as it comes out from 
the model calculation, before the support correction is implemented. For this reason it 
does not go to zero at x — 1. 

Figure 3: Proton OAM distributions in a modified IK model (see text) for the quarks, 
(a), and the gluons, (b). The full curves show the initial distributions at the hadronic 
scale of fj,Q = 0.23 GeV 2 , where around 40 % of the nucleon momentum is carried by the 
gluons; the dashed curves represent the LO-evolved distributions at Q 2 = 10 GeV 2 ; the 
long-dashed curves give the LO-evolved distributions at Q 2 = 1000 GeV 2 . 

Figure 4: Proton OAM distributions for the "D- Model" (see text), for the quarks, (a), 
and the gluons, (b). The full curves show the initial distributions at the hadronic scale of 
/2q = 0.23 GeV 2 , where around 40 % of the nucleon momentum is carried by the gluons; 
the dashed curves represent the LO-evolved distributions at Q 2 = 10 GeV 2 ; the long- 
dashed curves give the LO-evolved distributions at Q 2 = 1000 GeV 2 . 

Figure 5: The contributions to the proton spin sum rule, Eq. (|IJ), according to: (a) 
the modified IK scenario of Fig. 3; (b) the "D model" scenario of Fig. 4. The dashed 
curve shows |AS(Q 2 ), the long-dashed one Ag(Q 2 ), the dot-dashed curve is L q {Q 2 ), the 
dot-long-dashed curve gives L g {Q 2 ) and the full curve represents the sum of the previous 
four terms, giving the spin sum rule (J = |). 
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